Bridging the Gap Between the Mode Coupling and the Random First Order 
Transition Theories of Structural Relaxation in Liquids 
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A unified treatment of structural relaxation in a deeply supercooled glassy liquid is developed 
which extends the existing mode coupling theory (MCT) by incorporating the effects of activated 
events by using the concepts from the random first order transition (RFOT) theory. We show 
how the decay of the dynamic structure factor is modified by localized activated events (called 
instantons) which lead to the spatial reorganization of molecules in the region where the instanton 
pops up. The instanton vertex added to the usual MCT depicts the probability and consequences 
of such an event which can be derived from the random first order transition theory. The vertex is 
proportional to exp(— A/s c ) where s c is the configurational entropy. Close to the glass transition 
temperature, T g , since s c is diminishing, the activated process slows beyond the time window and 
this eventually leads to an arrest of the structural relaxation as expected for glasses. The combined 
treatment describes the dynamic structure factor in deeply supercooled liquid fairly well, with a 
hopping dominated decay following the MCT plateau. 
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Inelastic neutron scattering has given detailed struc- 
tural information on how transport in liquids changes 
as they are cooled in the glassy state P, Many of 
the structural details are well described by mode cou- 
pling theory (MCT) but that approach does not de- 
scribe well the dynamics of the deeply supercooled state 
H H H H 0, II The random first order transition 
(RFOT) theory of glass explains this lacuna of MCT 
[9L floL flll | . In this paper we show how RFOT theory can 
bridge the gap between the onset of glassy dynamics at 
high temperatures and the low temperature behavior and 
thus provide a unified view of the structural dynamics in 
deeply supercooled region. Over the last few decades, ex- 
perimental and theoretical studies on the phenomenolog- 
ical aspects of the glass transition have been the subject 
of a vast scientific literature. A glass forming liquid is 
characterized by its ability to circumvent crystallization 
on rapid cooling to a temperature well below its freezing 
temperature T m . In contrast to crystal formation what is 
most remarkable about solidification into an amorphous 
state is that the structure changes accompanying vitrifi- 
cation are either very small or unobservable. However, 
on passage through supercooled regime the dynamic and 
thermodynamic behavior of the liquid exhibits a number 
of anomalies as it eventually transforms to an amorphous 
solid, called a glass. 

The density-density correlation function at finite wave 
vector probed by neutron scattering, (pk{t) for a system 
approaching glass transition first decays via a fast mi- 
croscopic process which is followed by a plateau in the 
intermediate time. The dynamics while approaching the 
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plateau is non-exponential, power law in nature and is 
called f3 relaxation. The time scale of the (3 relaxation 
is predicted to show a power law divergence. The decay 
from the plateau known as the initial a relaxation is also 
given by a power law, known as von Schweidler law. The 
exponent in the (3 regime and the initial a relaxation are 
different in value but related to each other. The slow a 
relaxation which appears at long times is described by 
an stretched exponential or Kohlrausch- Williams- Watts 
(KWW) function Q. 

The thermodynamic and kinetic anomalies at T g are 
known experimentally to be correlated. One experimen- 
tally finds a sharp rise in the measured heat capacity of 
a liquid during heating which follows prior cooling at a 
constant rate[l2j, the cycle being well extended on cither 
sides of the glass transition region. The over shoot of the 
heat capacity is taken to be a thermodynamic signature 
of a glass to liquid transition. However, what lies at the 
heart of glass transition is the dramatic slowing down of 
molecular motion that occurs progressively on cooling af- 
ter the plateau in the inelastic structure factor appears. 
This is manifested by a phenomenal increase in shear vis- 
cosity r\ (also the characteristic structural relaxation time 
t) by several orders of magnitude for a relatively modest 
decrease of temperature (by few tens of degrees). The 
glass transition seen in the laboratory is in fact best de- 
scribed as a kinetic phenomenon which marks the falling 
out of the equilibrium due to the inability of molecular 
rearrangement within the experimental timescale as tem- 
perature T is lowered. A conventional definition of the 
transition temperature T g is the temperature at which 
the viscosity, r\ reaches a value 10 13 P. 

The mode coupling theory makes detailed predictions 
of the temperature dependence of the time evolution of 
the density-density correlation function Q. The emer- 
gence of a plateau in the neutron scattering is one of the 
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important predictions of the ideal mode coupling theory. 
At a quantitative level both the emergence of a plateau in 
the dynamic structure factor at high temperature and the 
exponential slowing of the dynamics at low temperature 
are also explained by the random first order transition 
(RFOT) theory 0, El In the high temperature limit, 
RFOT contains the essential elements of MCT as a stabil- 
ity limit. While idealized MCT can also correctly predict 
the power law dynamics in the f3 and the initial a relax- 
ation regime and the relation between the two exponents 
in these two regimes, it fails to explain the slow a relax- 
ation which appears at long times. MCT strictly applied, 
predicts a premature arrest of the decay of the density- 
density correlation function. The dynamical plateau of 
MCT was contemporaneously obtained by a variety of 
static approaches based on self-consistent phonon calcu- 
lations and density functional theory which form the ori- 
gin of the RFOT theory [H MCT would appear 
to predict a dynamic arrest of the liquid structural relax- 
ation, arising solely from a nonlinear feedback mechanism 
on the memory function. This feedback mechanism was 
first proposed by Geszti to describe the growth of 
viscosity in a previtrification region and involves a self- 
consistent calculation of the dynamical correlation func- 
tion and transport coefficients. In RFOT theory one sees 
the dynamic arrest is connected with the emergence of a 
mean field free energy landscape. 

According to the MCT, the essential driving force for 
glass transition is the slowing down of the density fluctu- 
ations near the wavenumber at q ~ q rn where the static 
structure factor is sharply peaked. Thus, prima facie the 
glass transition does not appear to be caused by a small 
wavenumber infra red singularity observable in the static 
structure factor, but rather it is a phenomenon where 
the intermediate wavenumbers are important, (although 
a diverging length associated with the four point correla- 
tions is expected in the mean field RFOT theory HHHI). 
There is a softening of the heat mode near q ~ q m . Kirk- 
patrick and Wolynes |l6( showed that this mode coupling 
transition was equivalent to the transition to an aperiodic 
crystal predicted by the quasi-static theories. Near the 
dynamical transition there is a general slowing down of 
all the dynamical quantities, and this slowing down is 
most effectively coupled to wavenumbers near q m . The 
glass transition temperature predicted by MCT always 
exceeds the laboratory T g . The explanation for this was 
pointed out by Kirkpatrick and Wolynes [l(| : Since there 
are exponentially many aperiodic structures any single 
one of them is unstable to transformation to the other. 
At T c the system is not frozen because activated events 
allow interconnection between these structures. Thus 
T c can be identified as a crossover temperature below 
which activated motions occur in which groups of par- 
ticles move in a cooperative way over local barriers in 
the free energy landscape. Computer simulation studies 
have amply demonstrated this change of transport mech- 
anism and have explicitly revealed the presence of such 
hopping motions in deeply supercooled liquids when the 



usual convective diffusive motions cease to exist [18|, [19| . 
Simulations also show that these hopping events relax 
the local stress in the system [20j . 

These properties of the activated events below T c can 
be represented by the random first order transition the- 
ory in a way superficially analogous to the theory of nu- 
cleation at random first order transition (although there 
are some differences). Calculation of fluctuations and 
transport in deeply supercooled or glassy medium are 
most easily made where the reference is an ideal glassy 
state and where the configurational entropy has vanished 
[H> HH E3l • The barriers for activated dynamics vanish as 
the temperature of the system approaches the dynamical 
transition temperature from below much like at a spin- 
odal. The barrier for activated events can be described as 
a nucleation phenomenon where the free energy change 
is a sum of the entropy change due to the formation of an 
entropy droplet and a surface term with a radius depen- 
dent surface energy 0, . A more explicit view of the 
similarities and the differences from nucleation appears 
in the recent work of Lubchenko and Wolynes and the el- 
egant analysis of Biroli and Bouchaud [ItJ ■ Deep in the 
supercooled regime, the nucleation free energy barrier is 
given by the following simple expression 

AF*/k B T = 32k B /s c (1) 

where s c is the configuration entropy of a single mov- 
ing bead i.e equivalent spherical particle. The coefficient 
32 arises from a specific microscopic calculation and re- 
flects the entropy cost of localizing a particle to vibrate 
in the cage, explicitly observed in the plateau of the dy- 
namic structure factor. If one approximates s c near the 
Kauzmann temperature by s c = Acp(T — Tk)/Tk, 
one then obtains the empirical Vogel-Fulcher-Tammann 
(VFT) law, with a relation between the liquid's fragility 
and the heat capacity jump, Acp. This expression agrees 
well with experimental results. However, Eq^] is an 
asymptotic form of the barrier, valid only near Tk and 
as will be discussed later, if we take into account the 
barrier softening effect near T c then the value of the 
barrier will be considerably smaller than that given by 

EqU] M- 

This effect is found to be stronger as we 
move away from Tk- Once the estimation of the bar- 
rier height is made the relaxation time r is calculated 
using, t/tq = exp(F* /ksT). The random first order 
phase transition theory also provides an explanation of 
the dominant source of non-exponentiality in relaxation 
dynamics in terms of distribution of the configuration en- 
tropy of the activated droplet. If one assumes a Gaussian 
distribution of the barrier, then the standard deviation 
of the barrier height at T g also turns out to be related to 
the specific heat discontinuity of the supercooled liquid. 
This correlation between nonexponentiality and AC P is 
observed in the laboratory [23|. 

In this article we use the concepts of mode coupling 
theory and random first order transition theory to in- 
clude both the convective diffusional motion (given by 
MCT) and the hopping motion (described by RFOT) in 
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the dynamics. Thus the present treatment also corre- 
lates the kinetic anomalies (like stretching of the relax- 
ation time) with the thermodynamic anomalies (jump in 
the heat capacity) primarily through the RFOT theory. 
Here we first discuss the effect of hopping alone on the 
structural density fluctuations. The rate of hopping and 
the size of the region involved in hopping are calculated 
from RFOT theory 10]. It is shown that hopping opens 
up a new channel for the structural relaxation. Next we 
present the equation of motion of the density-density cor- 
relation function where the diffusional and the hopping 
motions are considered as two parallel channels for the 
structural relaxation in the liquid. In the same spirit 
as in the case of the idealized MCT the new equation 
of motion which now has hopping motion in it, is cal- 
culated self-consistently with the longitudinal viscosity. 
The feedback mechanism which gave rise to the diver- 
gence of the viscosity in the idealized MCT Q is also 
present in our theory but due to the presence of hop- 
ping, the divergence is shifted to a lower temperature 
where the RFOT hopping rate now finally vanishes only 
as the configurational entropy goes to zero. In our pro- 
posed equation of motion the diffusional and the hopping 
motions are considered as parallel channels for the den- 
sity relaxation. In the first approximate scheme we write 
the total density correlation in the time plane and in the 
second approximate scheme it is written in the frequency 
plane. We show that not only both the approximations 
(although they appear to be quite different) lead to sim- 
ilar effects on the density-density correlation function, 
although, the predicted longtime decays are somewhat 
different. 



Before discussing our own theoretical scheme we would 
like to mention that Das and Mazenko |24| in their fluc- 
tuating nonlinear hydrodynamic theory and later Gotze 
and Sjogren 0, 0] in their extended MCT have ele- 
gantly demonstrated that in the relaxation kernel (mem- 
ory function), in addition to the density contribution, if 
a coupling to the current modes is also included, then 
below T c , the dynamic structure factor decays from the 
plateau value in the long time. Gotze and Sjogren 0, Q 
have shown that this modification eventually leads to a 
parallel decay channel of the dynamic structure factor 
which eliminates the ergodic to non-ergodic transition at 
T c . However, although the authors call this extra decay 
channel a hopping mode, there is no microscopic connec- 
tion made in those theories to the single particle or col- 
lective hopping or the system crossing over a free energy 
barrier in the free energy landscape. In those treatments, 
the strength and the timescale of this so called hopping 
mode was determined by fitting to experimental results, 
and an Arrhenius temperature dependence of the rate of 
"hopping" was assumed in this fitting. 



I. THEORETICAL SCHEME 

We propose a unified structural description cover- 
ing the whole temperature regime, from the high tem- 
perature microscopic dynamics described by collisional 
physics to the low temperature collective many particle 
activated dynamics. The mode coupling theory analyzes 
liquid state dynamics by starting from a gas-like mode 
of transport but including correlated collisions and back 
flow. MCT does a good job describing structural correla- 
tions observed by neutron scattering in the normal liquid 
regime. MCT can also describe the dynamics approach- 
ing the supercooled regime but fails below Tc , where ac- 
tivated processes also contribute substantially. A deeply 
supercooled liquid explores a very rugged underlying free 
energy surface. The dynamics within each well may still 
be described by the idealized MCT which does not take 
into account the possibility of hopping between wells. Ac- 
tivated events represent essential singularities in the noise 
strength and thus are not contained in the perturbative 
MCT. To include the effect of activated events on dynam- 
ical structure factor, we need to calculate two different 
quantities. First, we need to calculate the effect such 
an event (which we call "popping up of an instanton") 
would have on the dynamic structure factor. The impact 
of such an instanton depends both on its size, that is, 
how many molecules are involved in the activated event 
and on how far particles are displaced by an activated 
transition. In RFOT theory, these quantities are the dy- 
namical correlation length £ and the Lindemann length 
cIl ■ Second, we also need to know the probability of 
such an activated event occurring. We can estimate both 
this probability and the other quantities from the ran- 
dom first order transition theory using thermodynamic 
and structural information alone. 

A similar underlying physical picture to that presented 
here has been put forward by Fuchizaki and Kawasaki 
|25| using what they called dynamic density functional 
theory (DDFT) to go beyond the MCT. In this coarse- 
grained description one considers the time evolution of 
the density field on a lattice. The dynamics is supposed 
to be given by a master equation and the free energy 
functional is taken from Ramakrishna-Yussouff (RY) the- 
ory (2(|. As shown by Singh, Stoessel and Wolynes [Tij] . 
RY theory can produce multiple minima. Fuchizaki and 
Kawasaki have shown that transitions between such min- 
ima allow the structural relaxation of the coarse grained 
density field. They show these transitions give rise to 
a stretched exponential relaxation. In their calculation, 
the overlap between the time dependent density field and 
the reference one of the same lattice site, takes a constant 
value, close to unity, in the liquid but in the supercooled 
state, the overlap is less than one and decreases via dis- 
crete jumps which they have connected to activated hop- 
ping transition over free-energy barriers. They have fur- 
ther defined an order parameter which is the overlap of 
the density fields at different lattice sites. The distribu- 
tion of the order parameter, in the normal liquid regime 



4 



is sharply peaked around unity but at higher density the 
distribution is broad and peaked around a value less than 
unity. The broad distribution reflects the multi-minima 
aspects of the free-energy landscape. 

The scheme presented by FK was implemented numer- 
ically and does not give as explicit picture of the micro- 
scopic mechanism of hopping as the current treatment 
does. Moreover, due to the temporal coarse graining the 
initial decay of the dynamic structure factor, <fr(t) and 
the subsequent plateau are missed in their analysis. In 
addition, their analysis did not highlight the connection 
with the configurational entropy of the system. In con- 
trast the present unified theoretical treatment provides 
a 4>[t) which produces a correct behavior over the full 
time and temperature plane, and also provides an analyt- 
ical expression of the dynamic structure factor where the 
hopping alone is considered from a microscopic point of 
view. In the following section, we find 4>hop or the effect 
of hopping on the density relaxation. Subsequently we 
incorporate <fihop(t) in the total <f>{t) in a self-consistent 
way. 



II. EFFECT OF HOPPING ON DENSITY 
RELAXATION 

Here we describe the hopping in terms of the popping 
up of a single instanton which is an activated event whose 
probability is obtained from the random first order tran- 
sition theory. Let us say that an instanton pops up at 
a position R. We will take this region of influence to be 
spherical on the average. Near Ta, good arguments sug- 
gest the region actually is fractal or string-like [2^1 - But 
the effect will ultimately be spherically averaged anyway 
as we shall see below. The particles within the sphere 
of radius £ around the position R will be displaced by a 
Lindcmann length, di. Now the change in density p(r) 
due to a single instanton popping up can be written as, 

Sp new (r,t + St) = 6p(r,t) + 6((r- R) <0 



[ dt' [ dr' Sp(r,', t' 



x G(r,r',t-t')9((r' - R) < £) 



Sp(r,t) 



(2) 



Where 9((r — R) < £) provides the effect of the instanton 
felt at the position r, provided it is within the radius £. 
G(r, r't — t') is the Greens function which determines the 
effect of instanton in moving particles from position r to 
a new one r', typically Lindemann length away during 
time t and t + St. T)(R) determines the region where the 
effect of the instanton is felt. 

Since a given instanton pops up at a specific location 
the transformation is not translationally invariant, but 
when averaged over the droplet location and orientation 
we restore translational and rotational invariance, 



5p neW (r,t + St) = Sp(r,t) + ^ / ,/R.0((r \l !j 



f dt' [ dr'Sp{r',t') 

J JV(R) 

x G(r, r', t - t')9((r' - R) < £) 
-Sp(v,t) 



Sp(r,t) + ^ [ dt' [ dr'Sp(r',t' 
V J JV(R) 



•ir opyr '' ' 

IV(R) 

x G(r,r',t-t')n(r' -r) 



V 



Sp(r,t) (3) 



Where, 



fi(r-r') = /<2R0((r-R) < £)0((r' - R) < 
= ^(4e+|r-r'|)(2£-|r-r'|) 2 , 







= o 



for(r - r') < 2£ 
for(v - r') > 2£ (4) 



is the overlap volume of two spheres of radius £ centered 
at r and r', respectivelyand 



v = J dR0((r-R)<O = ^< 3 



(5) 



is the volume of the region participating in hopping. £ 
defines the re gion of hopping which is calculated from 
RFOT theory 10]. 

The RFOT £ is of the order of a molecular length, 'a' 
at T c but diverges as (T — Tk)~ 2 ' 3 as the Kauzmann 
temperature is approached. Near T g (where the timescale 
is of the order of 1 hour) £ is universally ~ 5 molecule 
spacings [Io| . 

Thus expanding Sp (r, t + St) about t we get, 



dt v 



dt' / dr' 5p(r',t') 

JV{R) 

: G(r,r',t-t')n(r' - r) 



v Sp(r,t) 



(0) 



Now lets say that the rate of a particle hopping is given 
by ph- Thus the probability of any particle hopping in 
the region of volume V during the time interval St is 
PhV x St/vp, where v p is the volume of a particle ( a 
particle will be a bead according to RFOT theory). If wc 
include this information in the above equation and take 
a Laplace transform we get, 
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s5p new (q, s) - Sp new (q,t = 0) 



Ph 



J dq 1 G(q 1 ,s)n(q-q 1 )5p(q, s) 



uo<5p(q, s) 



(7) 



Thus we obtain, 



</W(<7, s) 



s + (v Q - ClG(q, s))p h /v p 



(8) 



In calculation of the hopping rate we have considered that 
hopping is instantaneous. But simulations have shown 
[20j that once it occurs, hopping happens over a period 
of time. This would give a further frequency dependence 
to CI. The fluctuations of entropic driving force leads to 
a distribution of hopping rates ph 01 • To account for 
this we can consider that there are n-type of instanton 
with lifetimes (tj) and with occurrence probability p l h . 
In this case equation [5] can be averaged over different 
instanton lifetimes and rates. Here, for simplification we 
replace ph by the average hopping rate P which is given 
by P — —Hip l h . Since hopping rate is small (that is, hops 
are rare), we can assume that the duration or lifetime of 
a hop (instanton popping up) is mfinitesimally small. 
In the above equation, 



ttG(q,s) = y<fqiG(qi,*)n(q-qi) 



(9) 



In the above theory the length scale of the Greens func- 
tion is given by the Lindemann length di where as the 
length scale of f2 is given by £. The random first order 
transition theory gives us the estimate of £, that is the 
spatial extent of the instanton and also the rate at which 
the instantons pop up that is ph or more properly its 
distribution. The latter should be related to the waiting 
time distribution of hopping, as discussed by De Michele 
and Leporini ^ ■ 

Since £ >> di, G(q,s) has a much slower decay with 
q than fl(q — qi) which implies that the product decays 
slowly with q. The propagator G{q, s) should be just the 
short time part of the dynamic structure factor. 

Now we can write, 



<W0?, s) 



1 



(v -nG(q,s))P/v p 
1 



s + Khop(s) ' 



(10) 



As mentioned earlier, the length scale of the Greens 
function is determined by the Lindemann length di . If we 
neglect its frequency dependence then the Greens func- 
tion can be written as, G(q) — exp(—q 2 d 2 ). 



A. Calculation of the hopping rate from RFOT 
theory 

In this section we briefly discuss the calculation of the 
hopping rate P from the RFOT theory. As has been 
discussed before, the hopping rate of an instanton is con- 
nected to the free energy barrier, which in turn is con- 
nected to the configurational entropy of the system |9lll0| 
and an asymptotic form of that is given by EqQ] Note 
that, in glasses, we need the interface energy between 
two alternative amorphous packings and not between two 
thermodynamic phases, like a solid and a delocalized uni- 
form liquid. The alternative phases become more alike at 
T c . The intcrfacial energy (surface tension) that enters 
in the calculation of the instanton barrier is smaller at T c 
than at lower temperature. This effect has been discussed 
by Lubchenko and Wolynes who termed it as barrier soft- 
ening. A rigorous estimate of the barrier requires com- 
plex replica instanton calculations as recently outlined by 
Franz [23 and by Dzero, Schmalian and Wolynes |30| . 
An alternative method motivated by variational reason- 
ing proposed by Lubchenko and Wolynes is to write a 
simple interpolation formula for the free energy which 
allows the system to choose the smallest available bar- 
rier, and is given by [T(ij]. 



F(r) = 



T K T A 

t k + t a 



4?r ~ 
irr 3 Ts c 
3 



(11) 



where r k and Ta are the surface energy terms at Tk and 
Ta respectively. Ta gives an estimate of temperature at 
which a soft crossover to the onset of activated transport 
takes place or the temperature at which the barrier disap- 
pears when approached from below. The expression for 
the surface energy at Tk is, Tk = 4-7T a a 2 (r /a) 3 ^ 2 , where 
<7 = j(kBT/a 2 )ln((a/dL) 2 /ire) is the surface tension co- 
efficient and a is the molecular bead-size which refers to 
the size of the quasi-spherical parts of any a spherical 
molecule. The bead count 'b' for each molecule is also 
be calculated within RFOT theory. The expression for 
the surface energy at Ta is, Ta — 4-7T A fa 3 (r /a) 2 , where 

£{x = 4(^) 3/2 (l/ v ^+W^) ! where f = (1 - T/T A ) 
and A/, the excess free energy density computed from 
the density functional. Note that Eq^] tends to Tk 
as T — > Tk and similarly, to as T — ■> Ta- From 
the above expression of the barrier, the critical barrier 
free energy, F*(r), can be calculated as the maximum 
value of F(r). Microscopic calculation giving Ta can 
then be used with Eq^] to calculate the value of vis- 
cosity, r]jr\Q — exp(F* /fcgT), (where rjo is the viscosity 
at the melting temperature, T m and is obtained from 
experiments). Rather than carry out completely micro- 
scopic calculations, Lubchenko and Wolynes used this ap- 
proach to fit the experimental values of viscosity, where 
the values of the transition temperature Ta and Sfu 
(s c = Sfit(l —Tk/T)) are used as the fitting parameters. 
Tk is obtained from thermodynamics. For example, Tk 
is known to be 175K for Salol. 
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For our calculation of the hopping rate P, we use EalTTl 
with the fitted values of and Sfa from Lubchcnko 
and Wolynes 0. Thus P = ±exp(-AF* /k B T), where 
To is the timescale of the system at the melting point, 
T m which is taken from experimentally known value |3l| . 
The above expression would predict that rate of hopping 
becomes of the order of hours at T = T g and vanishes at 
T = T K 0. 

We also need an estimate of the size of the hopping 
region £ which has a temperature dependence and is 
equated to the r value where in EqEl F(r)=0. 



III. THE TOTAL DENSITY-DENSITY TIME 
CORRELATION FUNCTION 

In the last section we have evaluated the change in 
the dynamic structure factor due to hopping alone. In 
this section we incorporate this effect along with mode 
coupling feedback terms in the full density-density corre- 
lation function in a self-consistent manner. 

As mentioned earlier MCT describes well the dynamics 
for a moderately supercooled system. It also works well 
for a deeply supercooled system confined within a sin- 
gle free energy minimum. These features are preserved 
when the hopping term or instanton vertex is added. In 
the idealized MCT, the equation of motion for the dy- 
namic structure factor, (^-(-^(g, t) is first simplified by 
neglecting the coupling between different wave vectors 
and considering the contribution from the static and dy- 
namic quantities at a single wave number, q — q m , (q m is 
the wavenumber where the peak of the structure factor 
appears). The idealized MCT equation can be written as 
111, 

4 > MCT{t) + l^MCT^) + tia4>MCT(t) 

+ [ dt'^ MCT \t')^ MCT {t -t') = (12) 

Jo 

where §mct ^ s the normalized density time correlation 
function. This integro-differential equation has an un- 
usual structure but it can be solved numerically and 
many of its properties are now analytically understood. 
The frequency of the free oscillator can be approximated 
as, Q 2 = kBTq^ n /mS(q m ) and the damping constant, 7, 
which is the short time part of the memory kernel, is 
taken to be proportional to Slo- This is because random- 
izing collisions occur in Lennard- Jones liquids on nearly 
the same timescale as the vibrations. The fourth term on 
the left hand side has the form of a memory kernel and its 
strength is controlled by the dimensionless coupling con- 
stant which can be taken to be, A = (q m A 2 / 8n 2 p) S (q m ) , 
where A8{q — q m ) = S(q m ) — 1 @. It is useful to trans- 
form to the Laplace frequency plane, where the Laplace 
transform is defined as <p(s) = C[<p(t)]. Now the equation 
can be rewritten as, 



4>mct{s) = — njj— (13) 

s + s+nf'is) 

where the longitudinal viscosity rft is given by 

^( S )= 7 + 4A^£[^ CT 2 W] (14) 

The above equations are nonlinear in nature and when 
EqElis solved numerically or EasllHI and 1141 are solved 
self consistently for different values of A, a dynamical 
arrest of the liquid structural relaxation was predicted at 
T = T c 0,0,0- However, T c is always higher than T g . 
The origin for this inconsistency is the exclusion of the 
hopping motion in the pcrturbative MCT dynamics. 

Hopping and continuous motion are essentially two 
distinct channels for the structural relaxation in a su- 
percooled liquid, although they get coupled due to self- 
consistency, as described below. In the deeply super- 
cooled regime these two styles of motion are usually well 
separated in their timescales. Thus to include hopping 
in the dynamics of structural relaxation, as a first ap- 
proximation, valid in the deeply supercooled regime, the 
full intermediate scattering function can be written as a 
product of two functions |32( , 

<t>ii, *) - <f>MCT(q, t)(f>ho P (q, t) (15) 

where 4>Mcriq,t) is the mode coupling part and 
4>hop(q,t) is the contribution from the hopping motion, 
which is already discussed in the last subsection. 

We now derive the equation for 4>Mcriq,t) consistent 
with EqEl To begin with, note that the equation of 
motion <p Merit) should remain the same as 4> Merit) 
(EqO but due to the presence of hopping, the memory 
function will now be modified. The memory function, 
which is the longitudinal viscosity, is now determined by 
the full dynamic structure factor which is given by Ea ll5l 
Therefore, 4 1 Merit) is now solved self consistently with 
the full dynamic structure factor, <f>(t), and the equation 
of motion is written as, 

<P MCT{t) + l4> MCT it) + MCT it) 

+ 4\n 2 f dt'4> 2 {t')^MCTit -t') = o (16) 

Jo 

Note that for simplicity we have removed explicit q de- 
pendence in the above equation and all the quantities are 
calculated at q = q m . Although, the structure of EalT2l 
and Eg 1 161 are quite similar, the memory function in the 
former is calculated using only 4> Merit) whereas in the 
later it is calculated using the full intermediate scattering 
function, <p(t). Thus (j> Merit) is now dependent on the 
hopping motion. The equation in the Laplace frequency 
plane, for the modified 4>mct(s) remains the same as 
Eall3l except that the ideal longitudinal viscosity r\\ d is 
replaced by the following longitudinal viscosity given by, 
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% (s)=7 + 4AO2£[0 2 (t)] 



(17) 



We see at this order, the longitudinal viscosity is modified 
which in turn modifies the decay of <Pmct if) ■ 

Combining Ea 1 151 and Ea ll6l we can now write a equa- 
tion of motion for the full structure factor, <j>(t). The 
equation becomes simple when 4>hop{t) can be approxi- 
mated as an exponential with time scale \jKy lov . In the 
limit of small Khop the equation of motion for <p(t) can 
be written as, 

4>(t) + j<f>(t) + (K hopJ + n 2 )4>(t) 
+ A\n 2 f dt'4> hop (t')(f> 2 (t') 

Jo 

x 4>(t-t')+K hop ct>{t-t') =0(18) 

The effects of hopping motion in the full structural re- 
laxation can be incorporated directly in the frequency 
plane which is suggested by Eq""TO| and Eq[""""J| Exploit- 
ing the strict parallelism of hopping and convective mo- 
tion, these two equations can be combined to write an 
extended equation for the dynamic structure factor in 
the frequency plane as, 



1 



s = 



a + (wo - nG(q m ,s))P/v p + 
1 



s + Khop{s) + Kmct(s) 



(19) 



where 



Khop(s) = (v - Q,G{q- 
and 
Kmct(s) 



,s))P/v p 



fin 



(20) 



and r)i(s) is given by Ea ll7l 

Kmct(s) and Kh op (s) are the contributions to the rate 
of relaxation from MCT and hopping modes, respectively, 
and they act as two parallel channels for the decay of <p(i). 
At high temperature and low density Kmct » Khop, 
thus the decay of (f>(t) is determined primarily by the 
MCT term. But as we approach T c , the relaxation of 
<p(t) slows down which increases the longitudinal viscosity 
value (given by Ea ll4(l . This in turn causes Kmct(s) to 
decrease much faster than Kh op (s) at low frequency. Now 
if we neglect the frequency dependence of Kh op then the 
equation of motion for cj)(t) for small Khop value can be 
written as, 

4>{t) + 1 4>(t) + {K hopl + ttl)<t){t) 



AX® 2 



4>{t - t') + K hop q)(t - t') 



Note that as expected, while Eq[2]has similar struc- 
ture as Eq[""H| the memory function is different and the 
consequence of this will be discussed later. 

The second scheme (given by Ea ll9|) appears to be 
rather similar to extended mode coupling theory of Gotze 
and Sjogren^, Q. In the extended MCT the current 
contribution is incorporated in the memory kernel. The 
modified memory kernel leads to an extra decay channel 
of the dynamic structure factor thus successfully predict- 
ing the decay of <j){t) below T c . The rate of decay for the 
second channel, S, is called the rate due to the hopping 
mode. However, note that there is no direct connection 
between the expression of S in the Gotze-Sjogren anal- 
ysis and the hopping motion. Since the microscopic cal- 
culation of S can never predict a zero hopping rate they 
used 6 as a fitting parameter to explain the experimental 
results. 



IV. 



NUMERICAL CALCULATIONS AND 
RESULTS 



(21) 



For this paper we solved EqsEl HfillSI and |^ are 
solved numerically. Although in principle Eas 11912*01 can 

be solved self consistently with Ea ll7l in an iterative way, 
it requires a very large number of iterations to converge 
near the glass transition temperature. An alternative 
way is to numerically solve directly Eq[!"T] Both these 
calculations were done for a model system and for con- 
creteness were applied to a specific molecular system, Sa- 
lol. 



A. Results for the schematic model 

In the numerical calculations of the schematic model 
system, we took flo — 1, and 7 = f2o fixed. We varied A 
between .1-1, knowing that T c is reached for A = 1 and 
varied K hop from 0.01 to 0.0001. 

In figure 1 we plot both the idealized MCT result 
(given by Eall2|) and the modified full (j>(t) (given by 
Ea ll5|) dynamic structure factor against log(t). We see 
that although the ideal MCT saturates to a plateau value, 
the extended theory exhibits a hopping induced decay 
following the MCT plateau value - the duration of the 
plateau depends on the rate of hopping. 

In figure 2 we plot the total 4>{t) and 4>MCT(t) against 
time. The plot shows that hopping not only leads to the 
decay of the total (f>(t) (given by Ea ll5fl . but it also slows 
down the growth of the longitudinal viscosity thus facil- 
itating the structural relaxation even through the mode 
coupling channel. 

The long time behavior of ln(4>(t)) and ln(4>MCT(t)) 
are plotted against time in figure 3, for A = 1 and 
Kh op = -001. The semilog plot shows straight lines in 
the long time with different slopes for the MCT part and 
the total dynamic structure factor, thus indicating that 
both the functions are exponential in the long time but 




FIG. 1: The idealized MCT result (given by Eq[EJ and the 
modified full cj>(t) (given by EqEHJ have been plotted against 
log(t ps) for different values of Khop- All the plots are for 
A = 1 that is at T = T c , the mode coupling transition tem- 
perature. For different values of Khop the curve follows the 
idealized MCT result till it starts decaying from the plateau. 
The smaller the Khop value the longer is the plateau and the 
slower the long time decay. 




t (time) 

FIG. 2: The total </>(t) (given by Eq[TSJ, the modified 
<f>MCT{t) (given by Eq |16H and the idealized <f>MCT(t) (given 
by Ea ll2|l are plotted against time for A = 1 ( i.e T = T c ) 
and for Khop = .001. After the initial short time decay, both 
the total structural relaxation <f>(t) and the MCT part of the 
structural relaxation ^MCr(t) decays with time where as the 
idealized MCT theory predicts no decay of the structural re- 
laxation. This figure shows that the decay of <f>MCT{t) due to 
hopping induces continuous diffusion in the system. 



FIG. 3: ln(<j>(i)), ln(<pMCT(i)) and ln(<pMCT(t)) have been 
plotted against time where the parameters are the same as in 
figure 2. The plot shows that in the longtime both <f>(i) and 
4>MCT(t) decays exponentially with different timescales. Due 
to the presence of the explicit hopping term (see Eq |151 <j)(t) 
decays faster than cf>MCT(i)- 



with somewhat different time constants as evident in fig- 
ure 2. 

Next we compare the two different schemes at T = T c 
(that is, for A = 1) and for Khop — -001. The solu- 
tion of Eq^| and Eq|2] are plotted in figure 4 against 
log(t). Both the approximate schemes provide similar re- 
sults but the first scheme fEa ll8|) shows a faster decay 
implying that hopping has a stronger effect at this level. 
A comparison between Eq^| and Eq|2] shows that al- 
though the structures of the equations are similar, in the 
memory function of Ea ll8l there is an extra product with 
4>hop which eventually leads to a faster decay of the total 
function as seen in the figure. The advantage of the first 
approximate scheme is that we can separately calculate 
the MCT part of the dynamic structure factor and can 
then demonstrate explicitly how hopping induces a decay 
of 4>mct- This implies that hopping opens up continuous 
diffusion channels in the system. 

In the schematic calculation we have independently 
varied A and K) lop while keeping parameters Qo and 7 
constant. In real systems all of these are functions of 
density and temperature and needs to varied simultane- 
ously. Although the temperature dependence of Qq and 7 
are not too strong but A depends more strongly on T and 
Khop is expected to depend more strongly than exponen- 
tial on T . This means that as we increase the timescale 
of hopping (which is achieved by lowering the tempera- 
ture), the A value also increases. Thus for real systems, 
unlike as shown in figure 1, as l/Kh op gets stretched, the 
MCT plateau value also should increase somewhat. This 
is consistent with density functional theory of the aperi- 
odic solid. In the next subsection we apply our theory to 
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log 10 (t ps) 

FIG. 4: To compare the two different approximate schemes, 
give by Ea llSI (scheme 1) and by Eg 1211 (scheme 2) we plot 
4>{t) against log(t ps) for A = 1 (that is, at T = T c ) and 
for Khop = .001. The initial short time dynamics is identical 
for both the schemes but scheme 1 predicts a slightly faster 
longtime decay (see text for detailed discussion). 



TABLE I: In this table we present the value of the param- 
eters calculated. Qo is the frequency of the free oscillator, 7 
is the short time part of the memory function, acting as a 
damping constant. A provides an estimate of the strength of 
coupling and is a dimensionless quantity. £ is the radius of 
the region of hopping, ( v o~ Q G(q m )) gj veg tne e ff ec t Q f a s i n gi e 

hopping on the density and P is the hopping rate per particle. 
Although all the parameters vary with temperature, the tem- 
perature dependence of We find that S7o and 7 have a weak 
temperature dependence which implies that the short time 
dynamics remain unchanged. Although, A also shows a weak 
temperature dependence but a small variation in A leads to a 
substantial change in the dynamics. The hopping rate varies 
strongly with temperature. 



T (K) {lt(ps~ 2 ) jips- 1 ) A i/a P (ps- 1 ) 

270 1.353 1.163 0.94 1.89 L79 3.09 xlO -4 

256 1.268 1.126 1.01 2.32 3.47 3.2 xlO -7 

247 1.216 1.102 1.05 2.62 5.12 2.89 xl0~ 8 



a real system chosen as Salol and include these effects. 



B. Salol 

Salol (phenyl salicylate) is a fragile glass forming sub- 
stance. It is intra-molecularly hydrogen bonded, a van 
der Walls system which is considered as a model sub- 
stance for the study of glass transition and molecular 
mobility in the supercooled liquid. To perform the mi- 
croscopic calculations, we need to map the system into a 
Lennard-Jones system and for that we need the molecular 
diameter and the well depth. Although the shape of the 
molecule is not spherical we can calculate a hard sphere 
diameter a hs = 7.221 [33|. Note that in the RFOT the- 
ory the elementary particles in the system are not the 
molecules but beads and each molecule is made up of 
certain number of beads. We thus need to estimate the 
bead size a. Use of the RFOT theory (with softening) for 
the fit of the viscosity leads to the number of beads in a 
Salol molecule to be 6.29 [lTJ. Equating the volume of 
the molecule with the volume of 6.29 beads, we get the 
hard sphere bead size ahs = 3. 095 A which for the present 
calculation we equate to the Lennard-Jones bead size a. 
Next we need to estimate the temperature scaling which 
will be equivalent to the Lennard-Jones well depth e. To 
obtain the temperature scaling we use the fact that for 
reduced Lennard-Jones system the mode coupling transi- 
tion is known to take place at reduced density, p* = 0.95 
and reduced temperature, T* — 0.57. We are not aware 
of any reported density value of the Salol system, but a 
reduced temperature T* = 0.57 seems to be a reasonable 



value when compared with simulation results |ljj ■ For a 
Salol system it is known that the mode coupling transi- 
tion temperature, T c — 256K 34]. Thus the temperature 
is to be scaled by, e/k B = 449.122. 

Now that we know the molecular and the thermo- 
dynamic parameters of the Lennard-Jones system, we 
can calculate the various parameter values in the inte- 
gral equation that we need for a microscopic calculation, 
namely, for the MCT part, f2§ = k B Tq 2 n /mS(q m ), 7 and 
A = (q m A 2 /8ir 2 p)S(q m ) 1 where A5(q-q m ) = (S(q m )-l) 
|5j • Instead of microscopically calculating the short time 
part of the memory kernel, 7, we again take it to be pro- 
portional to Qo with the proportionality factor assumed 
to be unity (as our schematic model study). Gotze and 
Sjogren |g suggest that this proportionality factor even 
when varied between 1 and 100 does not effect the long 
time behavior. In the calculation of Q and A, S(q m ) is 
calculated for the above mentioned Lennard-Jones sys- 
tem. 

Next for the physical quantities in the hopping part of 
the calculation, we need to calculate £ and P. These are 
temperature dependent quantities. For this calculation 
we use Ea 1TT1 where T A = 333 and s fu = 2.65 10]. The 
values for all the calculated parameters are given in Table 
1. 

The plot for the Salol system, obtained by using the 
parameters reported in Table 1 is given in figure 5. The 
plot appears to be quite similar to figure 1, at and below 
the mode coupling transition temperature T c — 256K. 
Due to the presence of hopping the structural relaxation 
decays from the plateau value. From Table 1 we find that 
flo and 7 have a weak temperature dependence which is 
reflected in figure 5 in the near invariance of the short 
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log 10 (t ps) 

FIG. 5: The calculated structural relaxation for Salol at three 
state points which also includes the mode coupling transition 
temperature T c — 256K. Eq |18l is solved with the param- 
eters given in table 1 to obtain (j>(t). The solid lines A, B, 
and C correspond to <f>(t) at temperatures 270K, 256K and 
247K respectively. The dot-dashed lines are the plots at the 
same temperatures but without hopping (that is, the idealized 
MCT result). 



time dynamics at all the three temperatures. On the 
other hand, although A shows a weak temperature de- 
pendence but a small variation in A leads to a substan- 
tial change in the dynamics. Above A = 1 the increase in 
its magnitude leads to an increase of the plateau value. 
The hopping rate varies strongly with temperature and 
has a stronger effect on the dynamics above A = 1 where 
at lower temperature the plateau is stretched to longer 
times. Below A = 1 (or T > T c ), where (f>(t) can de- 
cay completely via the diffusive channel (MCT part), the 
presence of hopping leads to a slightly faster decay of the 
structural relaxation. 



V. CONCLUDING REMARKS 

Experiments show that the structural relaxation in a 
supercooled liquid exhibits rich dynamics over many time 
and length scales. Often the physical origin of these dif- 
ferent dynamics are quite disparate. To explain the dy- 
namics of the liquid over the whole temperature regime 
under a single theoretical framework has been a chal- 
lenge. The random first order transition theory which 
at high temperature also contains the essential elements 
described in the perturbative mode coupling theory [l6j 
has been successful in explaining the liquid dynamics 
around the laboratory glass transition temperature T g . 
The mode coupling theory, on the other hand has been 
able to explain the high temperature dynamics extremely 
well. The theory is also known to be accurate for short 



times in the supercooled liquid below T c . Below T c on 
long time scales the idealized MCT theory fails because of 
the exclusion of the activated hopping motion. As RFOT 
theory makes clear that hopping motions correspond to 
instantons which gives rise to essential singularities at 
Tk- In this work we have proposed a unified structural 
description of the liquid, covering the whole temperature 
regime, from above T c to Tk- 

Towards this goal we first treated the effect of hopping 
on the density fluctuation. In the theory the effect of a 
single particle hopping is connected to the neighboring 
density fluctuations through a propagator. The typical 
rate of a particle hopping and the spatial extent of a 
hopping event are obtained from RFOT theory. After 
incorporating all the informations it is found that hop- 
ping acts as a channel for the decay of the structural 
relaxation (Ea llOfl . Distribution of hopping times are es- 
sential for the non-exponentiality of the a relaxation [Tl| 
and will be incorporated in future work. 

The total structural relaxation after incorporating the 
hopping motion was calculated using two different ap- 
proximate schemes. As in idealized MCT, the density 
relaxation is calculated self consistently with the longi- 
tudinal viscosity. In the first scheme (given by Ea ll6fl the 
structural correlations in the time plane are taken as a 
simple product of the MCT part and the hopping part. It 
is found that due to the inclusion of the hopping motion 
the arrest of the structural relaxation at T c , as predicted 
by the idealized MCT, disappears. The theory with thus 
modified total structural relaxation predicts a hopping 
dominated decay following the MCT plateau. In the first 
scheme the part of relaxation of <p(t) due to the diffusional 
motion (MCT part) is calculated separately (Ea lltjfl . but 
self-consistently with the full dynamic structure factor 
(which now includes the hopping motion). The theory 
predicts that hopping not only relaxes the total density 
correlation but also invokes a decay of the MCT part. 
This implies that hopping opens up continuous diffusion 
channels in the system. This is consistent with the co- 
existence of both kinds of motions observed in computer 
simulation 1 19|. 

It is interesting to consider the predictions of the 
present generalized theory regarding the temperature de- 
pendence of viscosity. The decay of the MCT part of the 
dynamic structure factor already predicts a slower rise of 
the viscosity with the temperature than the ideal MCT. 
The interesting point to note is that the hopping rate 
itself also decreases with the lowering of temperature. 
Thus, the dynamic structure factor now not only has a 
slow exponential decay in the long time, but the plateau 
also gets stretched to longer times. Both the slow expo- 
nential relaxation at long times and the stretching of the 
plateau at 'intermediate times' contribute to the increase 
of the viscosity. However, the resulting temperature de- 
pendence can be quite different from the prediction of the 
ideal MCT. In particular, viscosity is dominated by the 
power law part at low supercooling but at large super- 
cooling even the duration of the plateau and the power 
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law decay are determined mostly by the hopping rate. ported in parts from NSF (USA) and DST (India). 
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